Scaling of cosmological domain wall networks with junctions 
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We describe the results of the largest and most accurate three-dimensional field theory simulations 
of domain wall networks with junctions. We consider a previously introduced class of models which, 
in the limit of large number TV of coupled scalar fields, approaches the so-called 'ideal' model (in 
terms of its potential to lead to network frustration). We consider values of TV between TV = 2 
and TV = 20. In all cases we find compelling evidence for a gradual approach to scaling, with the 
quantitative scaling parameters having only a mild dependence on TV. These results strongly support 
our no-frustration conjecture. 



INTRODUCTION 

If our current understanding of particle physics and 
unification scenarios is correct, defect networks must nec- 
essarily have formed at phase transitions in the early uni- 
verse The type of defect that forms and its specific 
properties depend on the particular details of each sym- 
metry breaking, so a wide range of possibilities exist, with 
correspondingly different cosmological consequences Q . 

Domain walls are usually pathological 0], but it has 
been claimed 0] that if a domain wall network is frozen 
in comoving coordinates (or 'frustrates', as is often col- 
loquially put) then it can naturally explain the observa- 
tional evidence that points to a recent acceleration of the 
universe. This happens because the equation of state of 
a domain wall gas is given by w = p/p — —2/3 + v 2 , v be- 
ing the root-mean squared (RMS) velocity of the domain 
walls. Note that we require that w < — (1 + fi^/O^) < 
— 1/2 in order to accelerate the universe at the present 
time and consequently v needs to be small. In prac- 
tice the characteristic scale of the network, L, also needs 
to be tiny in order not to give rise to exceedingly large 
CMB fluctuations. The simplest domain wall models are 
known to reach a scaling regime (until they dominate the 
energy density of the universe) , as first pointed out in [j| 
and recently studied in detail in [1, 0, 0] ■ Nevertheless, 
it was thought that more complicated models, notably 
those having junctions, would eventually frustrate. 

In previous work [^, [1(3] we have studied the dynamics 
of domain wall networks with junctions, and investigated 
in detail energy, geometrical and topological constraints 
on the properties of domain wall networks. This led us 
to develop an 'ideal' class of models which includes, in 
the large TV limit, what we called the 'ideal' model (that 
is, the best candidate for frustration). However, a series 
of analytic and numerical arguments led us to a no frus- 



tration conjecture: even though one can build (purely by 
hand) special lattices that would be locally stable against 
small perturbations 2, 11 , [lj| , no such configurations 
are expected to ever emerge from any realistic cosmo- 
logical phase transition. Our high-resolution numerical 
simulations of the ideal and other models showed clear 
evidence of a gradual approach to scaling, which was sub- 
sequently confirmed in [131 ] . 

Still, in our previous work 0, [l(| we only considered 
numerical simulations in two spacetime dimensions. The 
present work is the third in this series of papers, and 
its goal is to eliminate this shortcoming. We report on 
the results of a series of massively parallel domain wall 
network numerical simulations of the 'ideal' class of mod- 
els in three spatial dimensions, which provide conclusive 
evidence for a gradual approach to scaling and hence 
strongly support our no frustration conjecture. We start 
by discussing various dynamical issues relevant for the 
evolution of domain wall networks and briefly explain 
the 'ideal' model. We then briefly describe our numeri- 
cal implementation. Our code is based on the algorithm 
of Press, Ryden & Spergel [f| , though with a number of 
crucial improvements that we will point out. Finally, we 
discuss our results and draw some conclusions. 



THE IDEAL DOMAIN WALL MODEL 



In [6| we introduced a phenomenological one-scale 
model for the evolution of domain wall networks that 
has been shown to provide a good approximation to the 
evolution of two key network parameters: the character- 
istic scale of the network, L, and the RMS velocity of the 
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domain walls, v. The evolution equations are 

(1) 

where iT is the Hubble parameter, c w is the energy loss 
efficiency, k w is the curvature parameter and we have de- 
fined a damping length scale, l/£d = 3/f + 1/^/, which 
includes both the effects of Hubble damping and parti- 
cle scattering. The characteristic scale of the network is 
defined as L = a/p where p is the average density in do- 
main walls and a is the wall mass per unit area. Note 
that if domain walls are an important contribution to the 
dark energy then p must be of the order of the critical 
density, p c , at the present time. 

If we ignore, for the moment, the effects of friction and 
the energy loss by the network (by making £f — > oo and 
c w = 0) it possible to show that a linear scaling solution 
is possible for a oc t a with a > 1/4. In the radiation 
era we obtain L = \/4/3k w t and v — l/y/3 while in the 
matter era we have L = y/3/2k w t and v = l/y/6. We see 
that in both eras we have L ~ k w t and relatively large 
velocities. If we require CMB temperature fluctuations 
generated by domain walls on scales of the order of Hub- 
ble radius to be smaller than 10~ 5 then one would need 
(LHf/ 2 < 1CT 5 or equivalently L < 1 Mpc (our conser- 
vative estimate in 0]). However, CMB observations im- 
ply that the fluctuations generated by the domain walls 
have to be smaller than 10 -5 down to much smaller scales 
(~ iJ _1 /100). This means that current constraints on L 
are expected to be roughly 2 orders of magnitude stronger 
(L < 10 kpc) which implies a very small curvature param- 
eter k w < 10~ 6 ). This clearly shows that the simplest 
domain wall scenario without junctions (with k w ~ 1) is 
ruled out as a dark energy scenario. It is easy to show 
that allowing for a non-zero c w leads to a larger L and 
consequently it does not help frustration On the 

other hand, including friction also does not help much, 
due to the limited amount of energy with which domain 
walls can interact conserving energy and momentum Q . 

This failure of the simplest domain wall scenario led us 
to consider more complex scenarios with junctions and in 
Q we investigated in detail energy, geometrical and topo- 
logical considerations that severely constrain the proper- 
ties of domain wall networks. In particular, in the con- 
text of 2D domain wall networks, we have shown using 
local energy considerations that two edge domains are al- 
ways unstable and that three, four and five edge domains 
will be unstable if only Y-type junctions occur in a given 
model. We have also demonstrated that increasing the 
average dimensionality of the junctions, (d), leads to a 
decrease of the average number of edges, (x), per domain 
(in particular if (d) > 6 then (x) < 3 and consequently no 
equilibrium configurations will ever form). Also, allow- 
ing for domain walls with different tensions contributes 
to increasing the instability since the walls with higher 



tension will tend to collapse thus increasing the dimen- 
sionality of the junctions which, in turn, will in general 
lead to the production of further unstable two edge do- 
mains. 

Another important aspect, not yet discussed in pre- 
vious papers, is related to the fact that the stability of 
a given domain depends on global considerations (those 
depending on the configuration of the surrounding do- 
mains) as well as local ones (those associated with the 
domain itself) and we expect the global ones to become 
more important as we increase the dimensionality of the 
junctions or consider specific domains with a large num- 
ber of edges. In particular it is possible to show that 
a domain with three edges only survives the local sta- 
bility analysis by very little in the case where there are 
only X type junctions (the potential energy after the col- 
lapse would be at most about 10% larger than before). 
In this case, we expect that, in general, non-local effects 
will make a three edge domain unstable. As a result in a 
model with only X type junctions the only possible stable 
configuration is the one in which all the domains have the 
same number of edges, 4, which never occurs in the con- 
text of realistic domain wall network simulations. Note 
that we are assuming the junctions to be free through- 
out the paper. Otherwise their energy-momentum con- 
tribution could not be neglected, spoiling the dark energy 
properties associated with a static domain wall network. 

The above energy, geometrical and topological consid- 
erations led us to propose a class of models with N scalar 
fields and N + 1 vacua with the property that all possi- 
ble domain walls have equal tensions [10(. The 'ideal' 
model is obtained in the limit N — > oo. For large N the 
collapse of a single domain will only very rarely lead to 
the fusion of two of the surrounding domains. This is 
clearly a very desirable feature from the point of view of 
frustration (see Q). Also, by requiring all domain walls 
to have equal tensions we avoid another potential source 
of instability. 

A specific realization of this class of models with N 
scalar fields and a scalar field potential with N + 1 vacua 
was given in fl(il | 

N+l N 

^E^-^o) 2 withr? = 5>«-^) 2 ' (2) 

3=1 »=1 

where p,. are the N+l coordinates of the vacua of the 
potential. We have chosen pi . to be the vertices of an 
(N+ l)-dimensional regular polyhedra, and fixed the dis- 
tance between the vacua to be equal to the parameter 
ro. Given that in this model all possible domain walls 
have equal tensions, only Y-type junctions will form. Al- 
though it is possible to have only X type junctions in a 
model with 4 minima, this is in general not the case for 
a larger number of minima if all possible vacuum con- 
figurations exist in a given simulation. In particular this 
means that it is not possible to construct a generalization 
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of the 'ideal' model in which all junctions are of the X 

typ e - _ 

In [10| | we performed 2D domain wall network simu- 
lations of the 'ideal' class of models with N = 4 and 
N = 7. These results are useful since they are much 
simpler (hence easier to understand), containing many 
important features which are also relevant in higher di- 
mensions. However, realistic 3D domain wall network 
simulations are required in order to test our no frus- 
tration conjecture. In the next section we shall briefly 
describe these simulations. 
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MASSIVELY PARALLEL SIMULATIONS 

We performed high-resolution field theory numerical 
simulations on the UK Computational Cosmology Con- 
sortium's COSMOS supercomputer using a modified ver- 
sion of the algorithm of Press, Ryden and Spergel 
These are described in detail in 0, Q. The PRS al- 
gorithm [B[ modifies the domain wall thickness in order 
to ensure a fixed comoving resolution. We measure the 
domain wall velocities using an algorithm analogous to 
that described in Q which removes the radiated energy 
from the walls which otherwise would contaminate the 
estimate of the velocities. This is clearly an important 
advantage over previous velocity estimations (see for ex- 
ample [5]). We defined the domain wall as the region 
where V(<j>) > aV m&K (V max being the maximum of the 
potential and < a < 1) and we estimated the velocities 
as 
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where a dot represents a derivative with respect to confor- 
mal time and 7 = (1 — v 2 ) 1 ! 2 . The comoving characteris- 
tic scale of the network was estimated as L c = L/a ~ 5/ f 
where / is the volume fraction of the box with domain 
walls (that is with V(<f>) > a U max following our previ- 
ous definition) and 5 is the thickness of a static domain 
wall. We verified that, for 0.2 < a < 0.6, our results 
are almost independent of the threshold a, as long as the 
domain wall is sufficiently resolved. We assume initial 
conditions where the scalar field at each point in the grid 
is associated with a randomly choosen minimum of the 
potential. 

The PRS algorithm was then parallelized with 
OpenMP directives, and optimized for the shared mem- 
ory architecture of COSMOS. The results to be discussed 
below come from series of matter era simulations of 128 3 
and 256 3 boxes, for all values of N between 2 and 20, 
and with 10 different runs for each box size and number 
of fields. In addition, a single 512 3 box was simulated 
for all N between 2 and 20. A more extensive series of 
simulations studying, among other things, the evolution 



FIG. 1: The asymptotic values of v* and L c /n for the ideal 
class of models with N ranging from 2 to 20 (dotted and solid 
lines respectivelly) . In Figs. [1] and [2] the error bars represent 
the standard deviation in an ensemble of 10 simulations. 
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in other backgrounds, will be described in a subsequent 
publication. 

memory is approximately 
/1024 2 MB, so for 512 3 boxes 
required (for 20 fields). An 
file can also be produced at 
specified timesteps which can then be used to gen- 
erate animations, an example of which is available at 
http : //www. damtp . cam. ac . uk/ cosmos/viz/movies/evo2J25620_m 
As a benchmarking example, a 512 3 simulation with 3 
fields (requiring about 14.5 Gb of memory) takes about 
14 seconds per step on 16 processors and only 5 seconds 
per step on 32 processors (as the memory ratio becomes 
favourable) without box output, and a complete run 
takes just 85 minutes. For larger runs the scalability is 
good if one keeps the memory smaller than 1 GB per 
processor. The largest simulation we have performed, a 
512 3 box with 20 scalar fields and box outputs at every 
timestep, took about 2 hours and 15 minutes on 128 
processors. 



RESULTS AND DISCUSSION 

In Fig. [T] we plot the asymptotic values of u* and 
L c /rj (dotted and solid lines respectivelly) for the ideal 
class of models with N ranging form 2 to 20 (we end our 
simulations when r\ becomes equal to the box comoving 
size). The error bars represent the standard deviation 
in an ensemble of 10 simulations. As we increase N the 
asymptotic value of L c /r) decreases which is expected as 
we get closer to the ideal model. It is also significant 
that the differences between the successive N results fol- 
ic/?? become increasingly smaller for large N which is a 
clear indication that the results obtained for N = 20 are 
already close to the N — > 00 results. On the other hand, 
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FIG. 2: The scaling exponents, A, for all TV's between 2 and 
20, for the 128 3 , 256 3 boxes (dotted and solid lines respec- 
tivelly). The *'s represent the exponent measured for single 
512 3 boxes. Note that frustration would correspond to A = 1. 



we do not find any significant dependence the velocities, 
f*, with N. 

In Fig. [2] we plot the scaling exponents, A, defined by 
L c /i] oc r/~ x , for all values of N between 2 and 20, for the 
128 3 , 256 3 and boxes (dotted and solid lines respectiv- 
elly). The error bars represent the standard deviation in 
an ensemble of ten simulations and the *'s represent the 
exponent measured for single 512 3 boxes. The scaling ex- 
ponents were computed using the results from the second 
quarter of the dynamical range of the simulations. We 
see that A is slightly greater than zero which indicates 
that there are small departures a the scaling solution. 
The fact that as we increase the box size, thus evolving 
the simulations for a longer dynamic range, A gets closer 
to zero is a clear indication that the networks are slowly 
aproaching a scaling solution. We have also performed 
3D simulations in a model with only X-type junctions 
and N = 3 (4 minima - see [8| for specific realizations) 
and found no significant improvements with respect to 
the N = 2 case (3 minima) with only Y-type junctions. 



CONCLUSIONS 

In this letter we presented the most compelling evi- 
dence to date that domain wall networks can not be the 
dark energy. Note that in order to be able to rule out the 
domain wall scenario for dark energy a very large and rich 
class of models has to be analysed in detail. This lead 
us to develop a model best suited for frustration (the 
'ideal' model). We have shown that even this model fails 
to produce a frustrated domain wall network. Current 
observational constraints using cosmic microwave back- 
ground and supernova data already strongly disfavour 
w = —2/3 as the equation of state of a single dark en- 
ergy component 14 1. However, we should bear in mind 



that these results are dependent on strong priors. Even 
if we take them for granted and accept that domain walls 
alone cannot be the dark energy, they could still make a 
significant contribution. Our results, however, seem to 
exclude even that rather more contrived possibility. 
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